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A measurement of the cross section for the inclusive production of isolated photons by the CDF 
experiment at the Fermilab Tevatron collider is presented. The measurement covers the pseudora- 
pidity region |»y 7 | < 1.0 and the transverse energy range E^, > 30 GeV and is based on 2.5 fb _1 
of integrated luminosity. The sample is almost a factor of seven larger than those used for recent 
published results and extends the E^ coverage by 100 GeV. The result agrees with next-to-leading 
order perturbative QCD calculations within uncertainties over the range 50 < E'^ < 400 GeV, 
though the energy spectrum in the data shows a steeper slope at lower Ej,. 

PACS numbers: 12.38.Qk, 13.85.Qk 
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The measurement of inclusive prompt-photon produc- 
tion constitutes a test of perturbative QCD (pQCD) with 
the potential to constrain parton distribution functions 
(PDF), while avoiding the complications of jet identifica- 
tion and energy measurements. The photon cross section 
is also sensitive to the presence of new physics at large 
photon transverse energy. In high-energy pp collisions, 
photons are mostly produced via quark-gluon Compton 
scattering or quark-anti-quark annihilation. In addition 
to these processes, photons can also be produced through 
the fragmentation of outgoing partons, though this con- 
tribution is reduced when the photon is required to be 
isolated from other particles in the final state. An isola- 
tion requirement is necessary to suppress the background 
from energetic tt° and i] mesons. 

Previous experiments have presented results on inclu- 
sive prompt photon production compared to next-to- 
leading order (NLO) pQCD predictions [E S IE Hi- In 
this Letter, we present a new measurement of the inclu- 
sive cross section based on 2.5 fb _1 of integrated lumi- 
nosity collected by the CDF II detector at the Tevatron. 
Photons are required to be isolated, have a pseudorapid- 
ity |?7 7 | < 1.0, and a transverse energy Ej, > 30 GeV 
This measurement extends the Ej. coverage by 100 GeV 
compared to previous measurements |4j , and presents re- 
duced systematic uncertainties as a result of using an 
improved background subtraction method based on the 
shape of the isolation distribution. 

The CDF II detector Q is a general-purpose particle 
detector located at the Tevatron collider. The charged- 
particle tracking system is immersed in a 1.4 T mag- 
netic field aligned coaxially with the beam line, and 
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provides tracking coverage in the pseudorapidity range 
1 77 1 < 2.0. A central preradiator detector surrounds the 
tracking system and samples the electromagnetic show- 
ers that begin in the material in front of it. This detec- 
tor consisted of multi-wire proportional chambers at the 
beginning of Run II, and was upgraded to scintillation 
tiles in 2004. Scintillator-based electromagnetic (EM) 
and hadronic (HAD) calorimeters arranged in projective 
towers of size Arj x Acj) = 0.1 x 0.26 provide a coverage of 
\r)\ < 3.6. The energy resolution of the CEM calorimeters 
for photons and electrons is ct/Et = 13.5%/a/Et © 2%, 
where © represents sum in quadrature. The central 
(|?7| < 1.1) electromagnetic strip chambers (CES) are 
multi-wire proportional chambers embedded inside the 
EM calorimeter and positioned at a depth corresponding 
to the expected maximum of the longitudinal shower pro- 
file (6 radiation lengths). Anode wires and cathode strips 
measure 4> and z respectively, providing a 2 mm position 
resolution in each direction for 50 GeV electrons. 

The data are collected using a three-level on-line event 
selection system (trigger) that selects events with at least 
one energy cluster consistent with a photon in the final 
state. A photon cluster consists of one to three con- 
secutive calorimeter towers in the r\ direction. Photons 
are collected with two trigger thresholds in Ej,, 25 GeV 
and 70 GeV. In order to reduce contamination from neu- 
tral meson decays, the low Ej, trigger requires photon 
clusters to be isolated. The isolation requirement at the 
second level uses simple box patterns of towers, while 
at the third level the extra energy inside a cone of ra- 
dius R = \J (A0) 2 + (A77) 2 = 0.4 around the cluster is 
required to be less than 10% of the energy of the clus- 
ter (E clu ). The low Ej, trigger also requires the lateral 
shower profile of the CES cluster to be consistent with 
that of electrons, as measured in test beam data. 

The event selection requires the primary vertex z po- 
sition to be within 60 cm of the center of the detector to 
maintain the projective nature of the calorimeter towers. 
In order to suppress beam-related backgrounds, cosmic 
rays, and calorimeter noise, as well as leptonic W de- 
cays JtJ, the missing transverse energy of the event has 
to be less than 80% of the transverse energy of the lead- 
ing photon candidate. This requirement reduces these 
backgrounds to less than one percent each, while pre- 
serving almost 99% of the photon signal. Photon can- 
didates are required to be matched to a photon clus- 
ter and to be fiducial to the CES detector. Photons 
with additional CES clusters are rejected, since contri- 
butions from neutral mesons occasionally produce mul- 
tiple clusters. At most, one low transverse momentum 
track {p^ k < lGeV/c + 0.5%E^/c) is allowed to point 
to the photon candidate. The fraction of the energy of 
the photon in the hadronic calorimeter has to be small 
(E HAD /E EM < 0.055 + 0.00045£ dn /GeV), and, for 
events collected using the low Et trigger, the CES shower 
associated with the photon candidate has to have a pro- 



■5 



file consistent with test beam electrons [8fl . 

The transverse energy of the photon is corrected to ac- 
count for non-uniformities in the calorimeter response, 
and calibrated using electrons from reconstructed Z 
bosons. Photon candidates are required to have Ej. > 30 
GeV and to be isolated in the calorimeter, Ej =qa — Ej\ — 
Elf° < 2 GeV, where E^ =0A is the transverse energy in 
a cone of R = 0.4 around the photon. The isolation re- 
quirement reduces the background from neutral mesons, 
but also suppresses the photon signal coming from parton 
fragmentation processes. 

While the selection criteria remove the bulk of the 
neutral meson background, substantial contamination re- 
mains, mainly corresponding to fluctuations in the frag- 
mentation of jets, leading to neutral mesons that carry 
most of the parton energy. To subtract these isolated ir° 
and r\ mesons, the isolation distribution from the data 
(without the Elf° < 2 GeV requirement) is fitted with 
signal and background templates formed using Monte 
Carlo simulation. The isolation distribution is sensitive 
to the differences between prompt photons and back- 
ground: photons produce a well-defined peak at low iso- 
lation while neutral mesons present a flatter shape. The 
signal template is obtained from a pythia 6.216 pho- 
ton Monte Carlo sample and the background template 
is constructed by selecting photons from meson decays 
in a QCD PYTHIA sample. In both cases the underlying 
event model of pythia is tuned to CDF jet data (tune 
A These events are passed through a GEANt[TiI| 

simulation of the detector and subjected to the same se- 
lection requirements as the data. As an example, Fig. [1] 
shows the isolation distribution in data compared to sig- 
nal and background templates for photons in the region 
70 < El < 80 GeV. The raw E^° is corrected for leak- 
age effects and pile-up contributions which occasionally 
yields a negative E^°, but if the negative E^° bins are 
not included in the fit, the photon fraction changes by less 
than 2%. At high £j (> 200 GeV) the signal template, 
as extracted from the Monte Carlo simulation, does not 
describe accurately the signal peak. This is attributed 
to deficiencies in the details of the shower simulation in 
the calorimeter and, to a lesser extent, to the model of 
the underlying event [13] . An independent correction 
is applied to modify the signal templates and improve 
the fitting process. At low Ej\ no correction is neces- 
sary while at very high EZ the template is shifted by 
—0.5 GeV and its width reduced by 50%. The correc- 
tion, which is not applied to the background templates, 
only changes the final results by few percent, and is ac- 
counted for in the study of systematic uncertainties (see 
below). As a result of the fitting procedure, the photon 
fraction (J 7 ) is extracted from the measured isolation dis- 
tributions. Figure [2] shows the values for T as a function 
of EZ, including the systematic uncertainties discussed 
below. 

The raw inclusive differential cross section as 



a function of EZ is defined as da / '(dEZ drf) — 
{N 1 T)/(AEZ Arf (.trig £), where iV 7 is the number of 
photon candidates in a given EZ bin, AEZ (Ary 7 ) is the 
size of the EZ {rp) bin, (trig the trigger efficiency, and 
C is the integrated luminosity. The trigger efficiency is 
approximately 100% in the kinematic region of the mea- 
surement. The measured cross section is corrected for 
acceptance, efficiency of the photon selection, and reso- 
lution effects back to the hadron level [HI using a bin- by- 
bin unfolding procedure and a sample of prompt-photon 
events simulated with pythia. To avoid any bias on the 
unfolding factors due to assumptions about the true EZ 
spectrum, the photon Monte Carlo sample is re- weighted 
to match the measured spectrum. The resulting unfold- 
ing factors vary between 0.638 ± 0.003 and 0.69 ± 0.01 
with little Ej, dependence. 

A detailed study of systematic uncertainties is carried 
out [IH • The largest contribution to the total uncertainty 
at high EZ is caused by the 1.5% uncertainty on the pho- 
ton absolute energy scale, due to a small energy depen- 
dence in the energy ratio of simulated and data electrons 
from the Z mass peak. This introduces an uncertainty on 
the measured cross section that varies between 6% and 
13% as EZ increases. At low EZ the dominant uncer- 
tainty source is the photon fraction. Different methods 
are considered to construct signal and background tem- 
plates and extract T: signal templates are defined using 
electrons from Z decays in data instead of using Monte 
Carlo simulated events; very simple templates (two bins 
in isolation) are considered, to remove the details of the 
isolation distribution in the fitting procedure; the photon 
signal is then extracted using background templates with 
and without the -Ej-dependent correction that is applied 
to the signal templates. In addition, a completely differ- 
ent method [8j based on the shower profile of the photon 
candidate in the CES detector and the number of conver- 
sions in the material in front of the preradiator detector 
is used to determine the neutral meson background. As a 
result, a conservative systematic uncertainty that varies 
between 13% at low EZ and 5% at high EZ (see Fig. [2]) is 
assigned to the measured photon signal, covering the re- 
sults from all the alternative methods. A 3% uncertainty 
on the measured cross section, approximately indepen- 
dent of EZ, reflects uncertainties on the determination 
of the photon acceptance, and a 5% uncertainty on the 
measured cross section at low EZ accounts for the un- 
certainty on the CES cut efficiency. An additional 10% 
uncertainty in the photon isolation energy introduces a 
1% uncertainty on the measurement. Performing the un- 
folding procedure using unweighted Monte Carlo samples 
resulted in a less than 1% effect on the measured cross 
section. The different sources of systematic uncertainty 
are added in quadrature. The total systematic uncer- 
tainty varies between 15% at low and very high EZ, and 
8% at intermediate EZ. Finally, an additional 6% un- 
certainty due to the measurement of the integrated lumi- 



6 



nosity is considered 14 1. 

The measured inclusive isolated prompt photon cross 
section as a function of Ej, is presented in Fig. [3] and 
Table fl] The data are compared to NLO pQCD predic- 
tions as determined by the jetphox [l5| program with 
CTEQ6.1M PDF [16[, normalization, factorization and 
fragmentation scales set to Ej, and photon isolation re- 
quirement as for the data. Variations of the scales by a 
factor of two change the prediction by 15% at low Ej 
and 8% at high Ej. The uncertainty on the predictions 
due to PDF varies between 4% at low Ej and 13% at 
high Ej,, as determined using the Hessian method [l7j ]. 
In addition, we have evaluated the theoretical prediction 
using the MRST04 [lj| PDF, and hnd it well inside the 
experimental and other theoretical uncertainties. 

The theoretical prediction includes an additional cor- 
rection factor, Chad(Ej), to account for the presence of 
non-pQCD contributions from the underlying event and 
fragmentation into hadrons, that tend to increase the en- 
ergy in the isolation cone. Chad is estimated, using Monte 
Carlo generated events, as the ratio between the nominal 
Ej, distribution at the hadron level and the one obtained 
after turning off both the interactions between proton 
and antiproton remnants and the string fragmentation 
in the Monte Carlo samples. Two different sets of tuned 
parameters in pythia (tune a and DW [id ]) are consid- 
ered, and the mean effect Chad = 0.91 ± 0.03 , observed 
to have little Ej dependence, is taken as the correction. 
The uncertainty on Chad covers the results obtained with 
the different pythia tunes. As expected, the correction 
reduces the predicted cross section, since the presence of 
underlying event activity results in photons failing the 
isolation requirement. 

A difference between data and the NLO pQCD pre- 
dictions is observed for Ej < 50 GeV. Discrepancies 
were also observed at low px in previous measurements 
at collider and fixed target experiments [H, 0, For 
Ej, > 50 GeV, good agreement is observed, and a global 
X 2 test in this region, including correlations between sys- 
tematic uncertainties across Ej bins, finds a probability 
of 21%. 

In conclusion, the measurement of the cross section 
for the inclusive production of isolated prompt photons 
with 30 < E T < 400 GeV and |r?T| < 1.0 in pp col- 
lisions at y/s — 1.96 TeV, using data corresponding to 
an integrated luminosity of 2.5 fb , has been presented. 
The sample is almost a factor of seven larger than those 
used for recent published results 0| and extends the Ej 
coverage by 100 GeV. A new method to determine the 
photon fraction based on the shape of the isolation distri- 
bution was implemented for the first time at CDF, which 
resulted in smaller uncertainties compared to previous 
results [lj . The measured cross section agrees with NLO 
pQCD predictions within uncertainties except at low Ej 
where the data have a steeper slope. 
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FIG. 1: Measured isolation distribution for photons with 
70 < < 80 GeV. A x 2 fit to the data (full points) with pho- 
ton (dark histogram) and neutral mesons (light histogram) 
templates is used to extract the photon fraction. The result 
of the fit is shown as a full line with the associated uncertainty 
from the fitting procedure. 
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Elj,. The systematic uncertainty band is discussed in the text. 
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section for photons in the pseudorapidity region |?7 7 | < 1.0 
and 30 < E^, < 400 GeV. The uncertainties in the central col- 
umn are statistical. The additional 6% luminosity uncertainty 
is not included in the table. A parton-to-hadron correction 
(Chad = 0.91 ± 0.03) is applied to the pQCD predictions. 
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the theoretical predictions. 
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